Partitioned time discretization for atmosphere-ocean interaction

Numerical algorithms are proposed, analyzed and tested for improved efficiency and reliabil-ity of the dynamic core of climate codes. The commonly used rigid lid hypothesis is assumed,which allows instantaneous response of the interface to changes in mass. Additionally, mois-ture transport is ignored, resulting in a static interface. A central algorithmic feature is thenumerical decoupling of the atmosphere and ocean calculations by a semi-implicit treatmentof the interface data, i.e. partitioned time stepping. Algorithms are developed for simpli-fied continuum models retaining the key mathematical structure of the atmosphere-oceanequations. The work begins by studying linear parameterization of momentum flux in terms of windshear, coupling the equations. Partitioned variants of backward-Euler are developed allowinglarge time steps. Higher order accuracy is achieved by deferred correction. Adaptations aredeveloped for nonlinear coupling. Most notably an application of geometric averaging isused to retain unconditional stability. This algorithm is extended to allow different size timesteps for the subcalculations. Full numerical analyses are performed and computationalexperiments are provided. Next, heat convection is added including a nonlinear parameterization of heat flux interms of wind shear and temperature. A partitioned algorithm is developed for the atmo-sphere and ocean coupled velocity-temperature system that retains unconditional stability.Furthermore, uncertainty quantification is performed in this case due to the importance ofreliably calculating heat transport phenomena in climate modeling. Noise is introduced in two coupling parameters with an important role in stability. Numerical tests investigate thevariance in temperature, velocity and average surface temperature. Partitioned methods are highly efficient for linearly coupled 2 fluid problems. Exten-sions of these methods for nonlinear coupling where the interface data is processed properlybefore passing yield highly efficient algorithms. One reason is due to their strong stabilityproperties. Convergence also holds under time step restrictions not dependent on mesh size.It is observed that two-way coupling (requiring knowledge of both atmosphere and oceanvelocities on the interface) generates less uncertainty in the calculation of average surfacetemperature compared to one-way models (only requiring knowledge of the wind velocity)

Bayesian multiscale deconvolution applied to gamma-ray spectroscopy

A common task in gamma-ray astronomy is to extract spectral information, such as model constraints and incident photon spectrum estimates, given the measured energy deposited in a detector and the detector response. This is the classic problem of spectral “deconvolution” or spectral inversion. The methods of forward folding (i.e., parameter fitting) and maximum entropy “deconvolution” (i.e., estimating independent input photon rates for each individual energy bin) have been used successfully for gamma-ray solar flares (e.g., Rank, 1997; Share and Murphy, 1995). These methods have worked well under certain conditions but there are situations were they don’t apply. These are: 1) when no reasonable model (e.g., fewer parameters than data bins) is yet known, for forward folding; 2) when one expects a mixture of broad and narrow features (e.g., solar flares), for the maximum entropy method; and 3) low count rates and low signal-to-noise, for both. Low count rates are a problem because these methods (as they have been implemented) assume Gaussian statistics but Poisson are applicable. Background subtraction techniques often lead to negative count rates. For Poisson data the Maximum Likelihood Estimator (MLE) with a Poisson likelihood is appropriate. Without a regularization term, trying to estimate the “true” individual input photon rates per bin can be an ill-posed problem, even without including both broad and narrow features in the spectrum (i.e., amultiscale approach). One way to implement this regularization is through the use of a suitable Bayesian prior. Nowak and Kolaczyk (1999) have developed a fast, robust, technique using a Bayesian multiscale framework that addresses these problems with added algorithmic advantages. We outline this new approach and demonstrate its use with time resolved solar flare gamma-ray spectroscopy

Energetic proton spectra in the 11 June 1991 solar flare

We have studied a subset of the 11 June 1991 solar flare γ-ray data that we believe arise from soft proton or ion spectra. Using data from the COMPTEL instrument on the Compton Observatory we discuss the gamma-ray intensities at 2.223 MeV, 4–7 MeV, and 8–30 MeV in terms of the parent proton spectrum responsible for the emission

Teaching the Blind to Find Their Way by Playing Video Games

Computer based video games are receiving great interest as a means to learn and acquire new skills. As a novel approach to teaching navigation skills in the blind, we have developed Audio-based Environment Simulator (AbES); a virtual reality environment set within the context of a video game metaphor. Despite the fact that participants were naïve to the overall purpose of the software, we found that early blind users were able to acquire relevant information regarding the spatial layout of a previously unfamiliar building using audio based cues alone. This was confirmed by a series of behavioral performance tests designed to assess the transfer of acquired spatial information to a large-scale, real-world indoor navigation task. Furthermore, learning the spatial layout through a goal directed gaming strategy allowed for the mental manipulation of spatial information as evidenced by enhanced navigation performance when compared to an explicit route learning strategy. We conclude that the immersive and highly interactive nature of the software greatly engages the blind user to actively explore the virtual environment. This in turn generates an accurate sense of a large-scale three-dimensional space and facilitates the learning and transfer of navigation skills to the physical world

Energetic proton spectra in the 11 June 1991 solar flare

The June 11, 1991 gamma-ray flare seen by the Compton Gamma-ray Observatory (CGRO) displays several features that make it a dynamic and rich event. It is a member of a class of long duration gamma-ray events with both 2.223 MeV and greater than 8 MeV emission for hours after the impulsive phase. It also contains an inter-phase between the impulsive and extended phases that presents a challenge to the standard gamma-ray line (GRL) flare picture. This phase has strong 2.223 MeV emission and relatively weak 4.44 MeV emission indicative of a very hard parent proton spectrum. However, this would indicate emission greater than 8 MeV, which is absent from this period. We present the application of new spectroscopy techniques to this phase of the flare in order to present a reasonable explanation for this seemly inconsistent picture

Gamma ray measurements of the 1991 November 15 solar flare

The 1991 November 15 X1.5 flare was a well observed solar event. Comprehensive data from ground-based observatories and spacecraft provide the basis for a contextual interpretation of gamma-ray spectra from the Compton Gamma Ray Observatory (CGRO). In particular, spectral, spatial, and temporal data at several energies are necessary to understand the particle dynamics and the acceleration mechanism(s) within this flare. X-ray images, radio, Ca XIX data and magnetograms provide morphological information on the acceleration region [4,5], while gamma-ray spectral data provide information on the parent ion spectrum. Furthermore, time profiles in hard X-rays and gamma-rays provide valuable information on temporal characteristics of the energetic particles. We report the results of our analysis of the evolution of this flare as a function of energy (∼25 keV–2.5 MeV) and time. These results, together with other high energy data (e.g. from experiments on Yohkoh, Ulysses, and PVO) may assist in identifying and understanding the acceleration mechanism(s) taking place in this event

COMPTEL measurements of MeV gamma-ray burst spectra

We present results from the on-going spectral analysis of gamma-ray bursts measured by the COMPTEL instrument in its main Compton “Telescope” observing mode (0.75–30 MeV). Thus far, 18 bursts have been analyzed from three years (April 1991–April 1994) of observations. The time-averaged spectra of these events above 1 MeV are all consistent with a simple power law model with spectral index in the range 1.5–3.5. Exponential, thermal bremsstrahlung and thermal synchrotron models are statistically inconsistent with the burst sample, although they can adequately describe some of the individual burst spectra. We find good agreement between burst spectra measured simultaneously by BATSE, COMPTEL and EGRET, which typically show a spectral transition or “break” in the BATSE energy range around a few hundred keV followed by simple power law emission extending to hundreds of MeV. However, the temporal relation between MeV and GeV (e.g., as measured by EGRET) burst emission is still unclear. Measurement of rapid variability at MeV energies in the stronger bursts provides evidence that either the sources are nearby (within the Galaxy) or the gamma-ray emission is relativistically beamed

X- and gamma-ray observations of the 15 November 1991 Solar Flare

This work expands the current understanding of the 15 November 1991 Solar Flare. The flare was a well observed event in radio to gamma-rays and is the first flare to be extensively studied with the benefit of detailed soft and hard X-ray images. In this work, we add data from all four instruments on the Compton Gamma Ray Observatory. Using these data we determined that the accelerated electron spectrum above 170 keV is best fit with a power law with a spectral index of −4.6, while the accelerated proton spectrum above 0.6 MeV is fit with a power law of spectral index −4.5. From this we computed lower limits for the energy content of these particles of∼1023 ergs (electrons) and ∼1027 ergs (ions above 0.6 MeV). These particles do not have enough energy to produce the white-light emission observed from this event. We computed a time constant of 26+20−15 s for the 2.223 MeV neutron capture line, which is consistent at the 2σ level with the lowest values of ∼70 s found for other flares. The mechanism for this short capture time may be better understood after analyses of high energy EGRET data that show potential evidence for pion emission near ∼100 MeV

COMPTEL gamma-ray observations of the C4 solar flare on 20 January 2000

The “Pre-SMM” (Vestrand and Miller 1998) picture of gamma-ray line (GRL) flares was that they are relatively rare events. This picture was quickly put in question with the launch of the Solar Maximum Mission (SMM). Over 100 GRL flares were seen with sizes ranging from very large GOES class events (X12) down to moderately small events (M2). It was argued by some (Bai 1986) that this was still consistent with the idea that GRL events are rare. Others, however, argued the opposite (Vestrand 1988; Cliver, Crosby and Dennis 1994), stating that the lower end of this distribution was just a function of SMM’s sensitivity. They stated that the launch of the Compton Gamma-ray Observatory (CGRO) would in fact continue this distribution to show even smaller GRL flares. In response to a BACODINE cosmic gamma-ray burst alert, COMPtonTELescope on the CGRO recorded gamma rays above 1 MeV from the C4 flare at 0221 UT 20 January 2000. This event, though at the limits of COMPTEL’s sensitivity, clearly shows a nuclear line excess above the continuum. Using new spectroscopy techniques we were able to resolve individual lines. This has allowed us to make a basic comparison of this event with the GRL flare distribution from SMM and also compare this flare with a well-observed large GRL flare seen by OSSE

The angular distribution of COMPTEL Gamma-Ray bursts

The superior burst location capability of the COMPTEL instrument aboard the Compton Gamma-Ray Observatory allows us to study the small-scale angular distribution of burst sources with good sensitivity even though the number of burst detections is small. We accumulate four years (April 1991–April 1995) of observations to form a catalog of 27 burst locations whose mean 1σ uncertainty is ∼1°. We find that the COMPTEL bursts are consistent with an isotropic distribution of sources, yet the spatial coincidence of two of the bursts within COMPTEL’s angular resolution indicates the possibility of repetition. This possibility is studied using the two-point angular correlation function and the nearest neighbor statistic. Model dependent upper limits on the fraction of repeating sources are derived